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ABSTRACT 

Aims. Within the framework of the COSMOGRAIL collaboration we present 7- and 8.5-year-long light curves and time-delay esti- 
mates for two gravitationally lensed quasars: SDSS J1206+4332 and HS 2209-1-1914. 

Methods. We monitored these doubly lensed quasars in the R-band using four telescopes: the Mercator, Maidanak, Himalayan 
Chandra, and Euler Telescopes, together spanning a period of 7 to 8.5 observing seasons from mid-2004 to mid-2011. The pho- 
tometry of the quasar images was obtained through simultaneous deconvolution of these data. The time delays were determined from 
these resulting light curves using four very different techniques: a dispersion method, a spline fit, a regression difference technique, 
and a numerical model fit. This minimizes the bias that might be introduced by the use of a single method. 

Results. The time delay for SDSS J1206+4332 is Af^g = 1 1 1.3 ± 3 days with A leading B, confirming a previously published result 
within the error bars. For HS 2209H- 1914 we present a new time delay of AIba = 20.0 ± 5 days with B leading A. 
Conclusions. The combination of data from up to four telescopes have led to well-sampled and nearly 9-season-long light curves, 
which were necessary to obtain these results, especially for the compact doubly lensed quasar HS 2209+1914. 
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1. Introduction 



|Refsdal| ( |1964| l was the first to realize that the determination 
of the time delay between different images of a gravitationally 
lensed source can be used to derive the Hubble constant Hq, the 
expansion rate of the Universe, if the mass distribution in the 
lensing galaxy is known. The other way around, if a value for 
Hq is assumed, the time delay between the lensed quasar images 
allows one to study the mass of the lensing galaxy. 

Even if this method presents the advantage of being inde- 
pendent of other calibrations, it has often been subject to con- 
troversies both on time-delay values (e.g. on QSO 0957H-561 
see [Press et al.|[r992a|b[ iPeireFaLlfTgoel [Kundic et al.||1997l l 
and on the Hubble constant Hq derived from these time delays. 



this value seeming lower than the values derived by other tech- 
niques ( Burud|2001j ). A decade ago, time delays were measured 
on relatively short light curves of only a few seasons (e.g. IBurud 



let al.|2002||Patnaik & Narasimha f200r;'Corbett et al.|1996 '), us- 
ing only one time-delay extracting method at a time. This could 
lead to biased results and an underestimation of the error bars as 



* Based on observations made with the 1.2-m Swiss Euler telescope 
(La Silla, Chile), the 1.5-m AZT-22 telescope (Maidanak Observatory, 
Uzbekistan), the 2.0-m HCT telescope (Hanle, India), and the 1.2-m 
Mercator Telescope. Mercator is operated on the island of La Palma by 
the Flemish Community, at the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias. 
** Light curves will be available at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsa rc.u-strasbg.fr/viz-l 
jbin/qcat?J/A-HA/???, and on http://www.cosmograil.org. 



pointed out by Eulaers & Magain (2011 1. Nowadays, the time- 
delay approach is considered to be competitive and highly com- 
plementary to other cosmological probes, and it gives similar 
resuks for Ho ( |Oguri|2007^ Linder|2011[|Suyu|2012[ ). 

COSMOGRAIL is a project aimed at measuring time de- 
lays in most known gravitationally lensed quasars. The use 
of 1-m to 2-m class semi-dedicated telescopes has allowed 
us to acquire light curves spanning up to ten years of data. 
In combination with a homogeneous way of analysing data 
as described in Tewes et al. (2012a| accepted) and 'Tewes] 
|et al.| ( |2012b[ submitted), the collaboration is contributing 
to a significant rise in reliable time-delay measurements of 
gravitationally lensed quasars. The first results included time 
delays for SDSS J1650-H4251 (Vuissoz et aL] [2007] ), and 
WFI J2033-4723 ([Vuissoz et al. i 2008) and more recently for 
HE 0435-1223 ([Courbin et al.llIOTll i and RX Jl 133-1231 
( [Tewes et al.[|2012b submitted). In this paper, we present 
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time delays for two doubly lensed quasars: SDSS J1206+4332 
(12''06'"29'.65,+43°32'17'.'60; 72000.0) and HS 2209+1914 
(22''11'«30'.30,+19°29'12'.'00;72000.0). 

SDSS J1206+4332 was discovered by |Oguri et JL] ( [2005) . 
They found a z - 1.789 quasar whose lensed images are sep- 
arated by 2'.'90. The redshift of the main lensing galaxy is es- 
timated at z = 0.748. One or two other nearby galaxies might 
contribute to the lensing potential. 

During the Hamburg Quasar Survey [Hagen et al. ( 1999| l 
identified the system HS 2209+1914 as a bright z, = 1.07 
quasar having a galaxy neighbour detectable on their Schmidt 
plates. Nothing else was published on this object until Chantry 
et al. ( 2010| l, who obtained precise astrometry for the two lensed 
quasar images and the lensing galaxy based on the deconvolution 
of Hubble Space Telescope images. 

In Section 2 we summarize the data acquisition and reduc- 
tion techniques, and briefly explain the methods used to deter- 
mine the time delays. In Section 3, we present the results for 
the doubly lensed quasar SDSS J 1206+4332, including photom- 
etry and time delays. In Section 4, we take a closer look at 
HS 2209+1914 for which we give a first time-delay estimate. 
The last section is a summary of our results. 



2. Data acquisition and reduction 

2.1. Data acquisition 

Both lensed quasars were monitored for more than seven years 
on different telescopes. The 1.2-m Mercator telescope, located 
at the Roque de los Muchachos Observatory on La Palma, 
Canary Islands (Spain), is equipped with the MEROPE cam- 
era having a field of view of 6f5 by 6f5 and a scale of 0'.'19 
per pixel. COSMOGRAIL observations by the Mercator tele- 
scope stopped in December 2008. The remotely operated 2m 
Himalayan Chandra Telescope (HCT) started systematic mon- 
itoring for the COSMOGRAIL collaboration in 2007 and has 
continued up to now. It is equipped with the Hanle Faint Object 
Spectrograph Camera (HFOSC), a 2048x4096 pixels CCD cam- 
era with a field of view of 10' by 10' and a scale of 0'.'296 per 
pixel. 

These two telescopes provided the bulk of the data, and were 
complemented by data from the 1.5-m telescope at Maidanak 
Observatory, Uzbekistan. Most of the data were obtained with 
the SITE CCD camera, a 800 x 2000 array camera with a scale 
of 0'.'266 per pixel and a field of view of 3!5 x 8!9. During two 
seasons, a second camera (SI) was used with an array of 4096 x 
4096 pixels and a field of view of 18! 1 by 18! 1 on the same scale 
of 0'.'266 per pixel. 

For the lensed system HS 2209+1914, some data were added 
from the 1.2-m Euler telescope, located at La Silla, Chile. It is 
equipped with a 2048 x 2048 CCD camera having a field of view 
of ir by ir and a scale of 0'.'344 per pixel. 

At every observing epoch, and as long as there were no 
technical or meteorological problems, five dithered images were 
taken in the 7?-band with an exposure time of six minutes each. 

2.2. Data reduction 

To analyse COSMOGRAIL data in a consistent and homoge- 
neous way, we use the semi-automated pipeline that was de- 
scribed in ITewes et al. ( 2012b) (submitted) to carry out the re- 
duction of the CCD images. We summarize the main points after 
the bias subtraction and flatfielding of the frames. 



Four stars are used to construct a point spread function (PSF) 
per frame. Photometry of the sources is obtained through simul- 
taneous deconvolution using the MCS algorithm (Magain et al. 
|1998[ ), where the main feature is the finite final resolution, which 
avoids artefacts. The deconvolved images have a pixel scale of 
half the telescope's pixel scale and have a Gaussian PSF with 
a two pixel full-width-at-half-maximum. The flux of the point 
sources is allowed to vary from one frame to another, unlike the 
lensing galaxy, which is part of the numerical background that 
is held fixed throughout the frames. The simultaneous deconvo- 
lution of all frames constrains the positions of the quasar images 
very well, even if not all data have been obtained under optimal 
seeing conditions. 

An estimation of the photometric shot noise ctn, calculated 
for each quasar image and frame, is given byFl 



O-N = A /. + 



N,k.+R^ 



(1) 



with ft the flux in the quasar image (expressed in number of 
photons), Nsky the sky level of the exposure, R the CCD read 
noise (both in photons per pixel), and S the PSF sharpness given 

by 



Y^^PSFijf 



(2) 



i.e. the sum over all pixels of the squared normalized fluxes of 
the PSF pixels. 

The data points in the light curves are the median photome- 
try per night obtained by the simultaneous deconvolution of the 
images. Their error bars are the maximum of two kinds of error 
bars: an empirical error bar and an intrinsic one. Indeed, the to- 
tal error on the night cannot be smaller than the one given by the 
spread of the measurements in that night, nor can it be smaller 
than the error bar given by the photon noise combined with the 
renormalization error, so taking the maximum of both is justi- 
fied. They are generally very close to one another. 

The empirical error reflecting the spread of the measure- 
ments during one night is estimated by the standard deviation 
of the measurements per night normalized by the number of im- 
ages per night. Since our photometry is a median value per night, 
the standard deviation is estimated via the median absolute de- 
viation (MAD) estimator (Hoagl liTet al.| |l983) corrected by a 
scale factor for normal distributions. The MAD is defined as the 
median of the absolute values of the residuals from the data's 
median: 



MAD = t7iedi {\Xi - t7iedj(Xj)\) 



(3) 



It is related to the standard deviation via cr - K * MAD. In case 



of a normal distribution, K - 



1.4862, with (/!"' being 



(^-'(0.75) 

the quantile function of the normal distribution. 

On the other hand, the intrinsic or theoretical error, calcu- 
lated as the quadratic sum of the median shotnoise per night and 
the median renormalization error per night, gives a theoretical 
estimate of the error on the combined data point. That neither 
the empirical nor the theoretical error systematically dominate 
is a proof of the coherence of both approaches, and it justifies 
the idea of taking the maximum of both. 



' Based upon Heyer, Biretta, 2005, WFPC2 Instrument Handbook, 
Version 9.1, Chapter 6 
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Finally, we put together the light curves of the same object 
coming from different telescopes by optimizing the magnitude 
shifts between the curves of the same image from different in- 
struments. 



2.3. Time-delay analysis 

To avoid a bias due to the method, we use four different tech- 
niques for our time-delay analysis. Three of them, a spline fit, a 
dispersion minimization, and a method based on the variability 
of regression differences, are described and tested extensively in 
|Tewes et"aL] ( |2012^ (submitted). The fourth method, the numer- 
ical model fit, was described and applied in Eulaers & Magain 
( 2011| l. We quickly summarize the main features of each of 
them: 



1. The dispersion-like technique considers the time delay to 
be the shift between the light curves that minimizes the dis- 
persion between these data. Microlensing is modelled by 
polynomials up to the third order, but no model light curve 
for the overall quasar variability is constructed here. 

2. The regression difference technique searches for the time 
delay that minimizes the variability of the differences be- 
tween the time-shifted continuous regressions of the light 
curves. Microlensing variability is not explicitly modelled 
here. 

3. The free knot spline technique models the intrinsic quasar 
variability, as well as the microlensing variations, by fitting 
splines simultaneously to the light curves. The curves are 
shifted in time to optimize this fit. 

4. The numerical model fit (NMF) constiucts a numerical 
model for the quasar variability, together with a linear mi- 
crolensing trend for a given time delay. The optimal time 
delay is the one that minimizes the difference between the 
data and this numerical model. 

These four techniques are based on very different principles, 
so we expect them to be sensitive to different sources of error 
By applying all of them to our light curves, we minimize the 
bias that might be due to the choice of the method. 

We use a Monte Carlo approach for the time-delay uncer- 
tainty calculations. For the first three methods, these are based 
on simulated light curves with known time delays based on the 
spline fit as explained in |Tewes it al. (2012a ) (submitted). For the 
last method, we use the numerical model light curve to which we 
add appropriate Gaussian noise. 



3. SDSS J1 206+4332 

3.1. Deconvolution results and light curves 

For SDSS J1206+4332 we use data coming from three tele- 
scopes. At Maidanak this lensed quasar was observed from 
January 2005 until July 2008. During the 2007 season, the sec- 
ond camera SI was used. Mercator data span the period March 
2005 till December 2008, and HCT has continued monitoring 
this lens since May 2007 until now. The combination of these 
data results in seven complete observing seasons from 2005 un- 
til 2011. 

The object is unvisible from August until the beginning of 
November, a period of more or less 100 days. The mean temporal 
sampling during the visibility window is about one point every 
fourth night for Mercator and Maidanak data, and one point per 
week for HCT. In the overlapping 2008 season this gives more 



than one point every second night. Table[T]gives an overview of 
the optical monitoring for this quasar. One epoch corresponds to 
one data point in the light curve, which is the median value for 
that night as explained in Section |2] S is the temporal sampling 
of the observations expressed as an approximative average in 
number of days excluding the seasonal gaps. 

The four stars, designated 1 to 4 in Figure [T] are used to 
model the PSF needed for the simultaneous deconvolution. If 
these stars did not satisfy our quality criteria (due to e.g. cosmic 
ray hits), a second set of four stars was used: star 1 in combina- 
tion with star 5 to star 7. Five stars, 1 to 5, were also used as flux 
calibrators. 

Figure l2] shows that the simultaneous deconvolution of 11 09 
R-band Mercator images clearly reveals the lens galaxy between 
the lensed quasar images, as well as two other galaxies close 
to the system, which have also been identified by Oguri et al. 
(2005). Table l2]presents the relative astrometry of the system as 
determined from this deconvolution. It shows a slightly different 
position for the B image and the galaxy Gl in comparison with 



the positions in Oguri et al. ( 2005 1 




Fig. 2: Result of the simultaneous deconvolution of 1109 R- 
band Mercator images of J 1206+4332. The pixel size is 0'.'095 
(i.e. half the pixel size of the Merope CCD camera at the 
Mercator telescope), and the resolution is 2 pixels full-width- 
half-maximum. In addition to the main lensing galaxy Gl be- 
tween the lensed images A and B we clearly identify two more 
galaxies: G2 to the north of the system and G3 to the east. North 
is up and east is left. 



Object 


X (arcsec) 


y (arcsec) 


A 


0.000 ± 0.004 


0.000 ± 0.004 


B 


0.073 ± 0.005 


2.970 ± 0.004 


Gl 


-0.411 ±0.042 


1.817 ±0.012 


G2 


1.482 ±0.006 


5.959 ±0.010 


G3 


-1.902 ±0.027 


2.353 ± 0.021 



Table 2: Asti'ometry for SDSS J1206+4332 relative to the 
brighter quasar image A. The positive directions for x and y are 
west and north respectively. 



E.Eulaers et al.: XIII. Time delays of SDSS 11206+4332 and HS 2209+1914 
Table 1: Overview of the optical monitoring for SDSS J1206+4332. 



Telescope 


Location 


Instrument 


Pixel scale 


Field 


Time span 


Epochs 


S 


Mercator 1.2-m 


La Palma, Canary Islands, Spain 


MEROPE 


a' 193 


6:5 X 6:5 


Mar 2005 - Dec 2008 


196 


-4d 


AZT-22 1.5-m 


Maidanak, Uzbekistan 


SITE 


a'266 


3:5x8:9 


Jan 2005 - July 2008 


57 


^Ad 


AZT-22 1.5-m 


Maidanak, Uzbekistan 


SI 


^'266 


i8:ix i8:i 


2007 


14 


^Ad 


HCT2m 


Hanle, India 


HFOSC 


a'296 


10' X 10' 


May 2007 -July 2011 


115 


^Id 
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Fig. 1: Part of the field around SDSS J1206+4332. It is a combination of 29 frames taken at the 1.2-m Mercator telescope with a 
seeing < V!2 and an ellipticity < 0.12. The stars used to model the PSF are labelled 1 to 4. In case one of these stars did not have 
sufficient quality, three other stars labelled 5 to 7 were used. North is up and east is left. 



We obtain an image separation of 2'.' 98 between the quasar 
images A and B, which is slightly higher than Oguri et aL] ( |2005| l, 
but a lower value for the flux ratio F^/Fb = 1.31, for which 



Oguri et al. (2005 1 found 1.48. Our value is the result of the 



deconvolution and has not been corrected for the time delay yet. 
We come back to this in the next section. 

The light curves that have been obtained as described in 
Section l2] are shown in Figure l3] The 2007 and 2008 seasons 
include data from all three telescopes and show that they match 
and complement each other well. We distinguish clear features 
of quasar variability over periods as short as two months, for 
example in the 2006 season of the B curve. Visually, the light 



curves are compatible with a time delay around - 100 days. 
Finally, we also observe an evolution of the flux ratio between 
the quasar images, especially in the 2010 and 201 1 seasons. 



3.2. Time delay 



Based on their lens modelhng, Oguri et al. ( |2005[ l predicted a 
time delay for SDSS J1206+4332 of AtAB = 92.6 or 104.4 days 
depending on the influence of the secondary galaxy G2. With 
only one observing season, Paraficz et al. (2009i found a time 



delay of AIab - i^l6^t days. 



'004 
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Fig. 3: Light curves of SDSS J1206+4332 combining data from Mercator, Maidanak, and HCT, and spanning 7 observing seasons. 
Photometry is obtained from simultaneous deconvolution, and the lo" error bars are calculated as described in Sectionl2| 



Our time-delay analysis is based on the light curves pre- 
sented in Figure [5] and it uses the four different methods briefly 
explained in Section l2] Since three of these techniques (disper- 
sion, regression difference and spline fit) rely on optimizations 
of a first guess of the time delay, we test the robustness of the 
result by running 200 times these techniques on the same data 
set, while uniformly and randomly varying the initial delay in 
a range of ±10 days around the first guess. Figure [4 shows that 
in the case of SDSS J 1206+4332 we have monomodal distribu- 
tions for the spline and regression techniques, and a bimodal dis- 
tribution for the dispersion technique, which is inherently more 
sensitive to the variation of this entry delay. The results of this 
test confirm that our method parameters are well adapted to the 
data, so that we can use the mean value of these distributions as 
our final best time delay. 

To obtain realistic error bars, we applied these three methods 
to 1000 synthetic light curves with known time delays. The er- 
ror bars for these three techniques consist of a systematic and a 
random part, which then are summed quadratically to obtain the 
total error. Figure l5] illustrates both contributions: we plot the 
mean measurement error as a function of the true delay. The er- 



ror bars are the standard deviation of these measurement errors. 
To be conservative, the maximum measurement error is our sys- 
tematic error, and the maximum standard deviation our random 
error, which are then summed quadratically. 

The numerical model fit (NMF) is completely independent 
of the three other methods in how it determines the time de- 
lay and its associated error bars. We add normally distributed 
random errors with the appropriate standard deviation to the nu- 
merical model light curve and redetermine the time delay. This 
procedure is repeated 1000 times. The mean value of the time- 
delay distribution that we obtain in that way is considered to be 
the final time delay, and its dispersion represents the Icr error 
bar Because the errors added to the model, hence the final error 
bars on the time delay, depend on the photometric error bars of 
our data, we checked that our residuals (data minus model) were 
compatible with a normalized Gaussian distribution. 

One of the parameters to be optimized in the NMF method 
is the flux ratio corrected by the time delay, for which we find 
Fa/Fb - 1.45 for the complete data set, which is very close to 
the |Oguriet all ( [20051 ) value of 1.48. 
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Fig. 4: Distribution of time delays for SDSS J 1206+4332 by 
randomly varying 200 times the initial time delay needed for 
three of our methods. Only the dispersion technique shows a sec- 
ondary peak. 
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Fig. 5: Error analysis of the simulated light curves for 
SDSS J1206+4332 using data from Mercator, HCT, and 
Maidanak. We plot binned mean measurement errors and their 
standard deviations against the true delays of the synthetic 
curves. The systematic error is largest for the dispersion tech- 
nique, and random errors dominate for the regression and spline 
methods. 



We present an overview of the time delays, along with 1- 
cr error bars, obtained by the different methods in Table l3] first 
for the Mercator and HCT data, two sets that really complement 
each other, and then for the whole data set coming from the three 
telescopes. Adding Maidanak data to the Mercator and HCT 
light curve adds more scatter, but seems to constrain the model 
light curves better, which is reflected in larger error bars for the 
dispersion method that does not involve a model, but smaller 
ones for the regression difference technique, the spline fit, and 
the numerical model fit. Image A leads image B by At^s days. 
Our results coming from four basically very different time-delay 
estimation techniques are very consistent, and confirm the exist- 
ing estimation from Paraficz et al. ( 2009 ), while improving its 
error bars for the spline fit and the NMF techniques used on the 
complete data set. 



That the error bars for the dispersion technique, the regres- 
sion difference, and the spline fit are systematically larger than 
for the NMF is inherent to the way they are calculated. For these 
first three methods, the error bars are the quadratic sum of two 
maximum contributions, as described before, whereas the NMF 
does not take the systematic error into account (generally under 
one day), and considers the standard deviation of the time delay 
to be its random error without taking a maximum over models 
with slightly different time delays. On top of that, the random 
error contribution of the first three methods is the result of red 
noise, whereas the NMF only takes white noise into account. 

Since the time delays resulting from the four methods all 
converge towards the same value, we only have to decide on the 
size of the error bars, which reflect the different sensitivities to 
various aspects that are present in the data and the method. Given 
the low spread in time-delay values between the methods, and 
given the possibility that the error bars from the dispersion meth- 
ods, the regression difference technique, and the spline fit have 
been overestimated, while those for the NMF might be underes- 
timated, we think that a final time delay for SDSS J1206+4332 
of AtAB - 111. 3 ±3 is a reasonable and comprehensive value, 
which can be applied in lens modelling. 



4. HS 2209+1914 

4.1. Deconvolution results and light curves 

The doubly lensed quasar HS 2209+1914 was observed by 
four telescopes involved in COSMOGRAIL: the Euler tele- 
scope started observations in May 2004 and was joined by 
the Mercator telescope, which continued the monitoring until 
December 2008. In May 2007, observations of this lens also 
started at the Maidanak Observatory during two seasons, and 
they were joined by the HCT the same year. This last telescope 
has continued observations until now. Even without the end of 
the 2012 season, we dispose of a total of 8.5 seasons of observa- 
tions from 2004 to mid-2012. 

The lensed quasar can be observed from May to December 
suffering from a non-visibility window of four months. The tem- 
poral sampling varies with the seasons and telescopes: combin- 
ing data from Mercator, HCT, and Maidanak yields the best sam- 
pling of more than one point every second night in 2007 and 
2008. This drops to one point every fortnight on average in 2010, 
2011, and 2012 when we only dispose of HCT data. Table |4] 
gives an overview of the optical monitoring for this lens. One 
epoch corresponds to one data point in the light curve, which is 
the median value for that night, as explained in Sectionl2] 

Figure l6] shows the four stars, labelled 1 to 4, used to model 
the PSF for the simultaneous deconvolution. Three of these stars, 
1 to 3, were also used as flux calibrators. 

In order not to degrade our photometry by some bad expo- 
sures, we eliminated all images with a seeing > 3" or an ellip- 
ticity of the point sources > 0.3, as well as images with too high 
a sky level or deviating normalization coefficient. Simultaneous 
deconvolution of all remaining 2242 images reveals a faint ring 
like structure (see Figure iTl) that is very similar to the decon 



volved Hubble Space Telescope images shown in [Chantry et al. 



(2010) . A secondary galaxy to the west of the system is also vis- 
ible. 

However, the main lensing galaxy cannot be resolved, given 
the compactness of the system. Both the image separation of 
only 1'.'02 and the flux ratio after deconvolution of (Fb/Fa) = 
0.82 agree with|Chantry et al.|(|20T0|, who find 1'.'04 and 0.79 ± 
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Method 


Mercator +HCT 


Mercator + HCT + Maidanak 


Dispersion 


11 1.52 ±4.95 (4.57, 1.89) 


113.65 ±6.79 (4.12, 5.40) 


Regression 


110.89 ± 10.04(9.55,3.08) 


109.73 ± 8.28 (7.29, 3.92) 


Spline 


111.18 ±5.79 (5.00, 2.92) 


111.31 ±3.93 (3.43, 1.92) 


NMF 


111.87 ±0.96 


113.80 ±0.90 



Table 3: Summary of time delays for SDSS J1206+4332 obtained by the four different methods. The indicated error is the total 
error. We mention between brackets the random and the systematic error for the first three methods as explained in the text. 

Table 4: Overview of the optical monitoring for HS 2209+1914. 



Telescope Location 

Eulerl.2-m ESO La Silla, Chile 

Mercator 1.2-m La Palma, Canary Islands, Spain 

AZT-22 1.5-m Maidanak Observatory, Uzbekistan 

AZT-22 1.5-m Maidanak Observatory, Uzbekistan 

HCT 2m Hanle, India 



Instrument Pixel scale 



Field 



Time span 



Epochs 



C2 


a' 344 


ll'x 11' 


May 2004 - Oct 2006 


21 


MEROPE 


a' 193 


6:5 X 6:5 


Aug 2004 - Dec 2008 


191 


SITE 


a'266 


3:5x8:9 


May 2008 - Aug 2008 


32 


SI 


a'266 


i8:ix i8:i 


May - Oct 2007 
+ Sept - Oct 2008 


58 


HFOSC 


a'296 


10' X 10' 


Sept 2007 -Aug 2012 


136 




Fig. 6: Part of the field of view around HS 2209+1914. It is a combination of 49 frames taken at the Mercator telescope with a 
seeing < I'.'l. The stars used to model the PSF are labelled 1 to 4. All but the last one are also used to calculate the normalization 
coefficient. North is up and East is left. 



0.027, respectively, even if our flux ratio has not been corrected 
for the time delay yet. 

Applying the data reduction described in Sectionl2J we show 
the complete light curve in Figure |8] Euler and Mercator data 
complement each other well in the 2004 and 2005 seasons. The 
2007 and 2008 seasons show a very good match between the 
data from the Mercator, Maidanak, and HCT telescopes. Even if 
we can observe brightness variations in the light curve up to 0.2 
mag over the years, clear features within a season are unfortu- 
nately absent, which complicates the time-delay measurement. 
However, the slight rise in magnitude at the end of 2009 in the B 
curve pins down the order, with the fluctuations in the B image 
preceding those observed in A. 



4.2. Time delay 

No theoretical estimations or previous measurements of the time 
delay are available for this system. Our analysis is based on the 
light curves shown in Figure [8j and it makes use of the four 
techniques of which we recalled the basic principles in Section 
I2] As we explained in more detail in the previous section on 
SDSS J1206+4332, we first test the sensivity of the results to 
the initial guess of a time-delay value necessary for the disper- 
sion, regression, and spline fit methods. In Figure l9] we see that 
the regression difference technique and, to a lesser extent, the 
dispersion method show a second solution. We take this into ac- 
count because our final best time delay is the weighted mean of 
these different values. 
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Fig. 8: Light curves of HS 2209+1914 combining data from Mercator, HCT, Euler, and Maidanak and spanning nearly 9 observing 
seasons. Photometry is obtained from simultaneous deconvolution, and the Icr error bars are calculated as described in Sectionl2| 



The error bars for the dispersion, regression difference, and 
spline fit techniques are calculated in the same way as in the 
previous section. The analysis of the random and systematic er- 
rors, based on 1000 synthetic curves with known time delays, 
and presented in Figure 10 show very small systematic errors in 
comparison to the large random errors. 

Applying the numerical model fit to our light curve consti- 
tutes an independent test for the time-delay results. To have re- 
alistic error bars on our time delay, we have to check, as previ- 
ously explained, whether the residuals from our fit are compat- 
ible with a normal Gaussian distribution with appropriate mean 
and sigma. This is not the case for HS 2209+1914, so we mul- 
tiply the photometric error bars of the data by the appropriate 
positive factor to obtain a correct distribution. As mentioned in 
Tewes et al. (2012aJ (submitted), a rescaling of these error bars 
does not directly influence the uncertainties on the time-delay 
measurement for the dispersion, regression difference, and spline 
fit techniques, but it does of course for the NMF method as the 
noise added to the model light curve directly depends on the pho- 
tometric error bars of the data points. Larger photometric errors 



lead to a larger uncertainty on the time delay obtained by the 
NMF 

According to the telescopes involved, the NMF method 
finds a flux ratio corrected by the time delay ranging between 
Fb/Fa = 0.82 (Maidanak) and Fb/Fa = 0.85 (Mercator) and 
of Fb/Fa = 0.84 when applied to the complete data set. These 
values are only slightly higher than the ratio obtained after de- 
convolution mentioned previously. 

We present an overview of our time-delay results from the 
four different methods on the complete data set in TablelS] Image 
B leads image A by AIba days. We notice that in spite of very 
large error bars, especially for the dispersion and regression dif- 
ference techniques, the different time-delay estimators do con- 
verge to a common time delay of AIba ~ 20 days with image B 
leading A. 

The large error bars for two of our methods can be par- 
tially explained if we take a closer look at the contributions of 
the data from different telescopes to the time delay. We applied 
the NMF to five subsets of data: 1) Mercator only, 2) Maidanak 
only, 3) HCT only, 4) Mercator + HCT, 5) Euler + Mercator + 
HCT. Comparing the results of a single run of the method on 
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Fig. 7: The numerical background common to the simultaneous 
deconvolution of 2242 R-band images of HS 2209+1914 reveals 
a faint ring-like structure that is very similar to the deconvolved 
Hubble Space Telescope images shown in 'Chantry et al. ( 2010| l. 
A secondary galaxy to the west of the system is also visible. 
However, the main lensing galaxy cannot be resolved given the 
compactness of the system. North is up and east is left. 
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Fig. 9: Distribution of time delays for HS 2209+1914 obtained 
by running three of the time-delay estimation techniques on 
the complete data set 200 times, while uniformly and randomly 
varying the initial delay in a range of ±10 days around our first 
guess for the time delay. 
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Fig. 10: Error analysis of the simulated light curves for 
HS 2209+1914 using all data from Euler, Mercator, Maidanak 
and HCT. We plot binned mean measurement errors and their 
standard deviations against the true delays of the synthetic 
curves. The random errors largely dominate the systematic ones, 
and are the smallest for the spline fit method. 



these subsets and a Monte Carlo approach of 1000 runs on per- 
turbed model curves, both presented in Table |6] where image 
B leads image A by Af^A days, reveals interesting information. 
Mercator and HCT data do not converge to the same time delay: 
Mercator data point to one around a; 23 days, whereas HCT data 
allow a possible delay around x 33 days. With Maidanak data, 
which cover the two overlapping seasons between Mercator and 
HCT exactly, we find both delays. With bare Maidanak data, the 
lower Mercator value is preferred, but once we perturb the model 
curve, we find a bimodal distribution with peaks around s; 22 
and a: 34 days. The more data we add to the light curve, the 
more the lower delay dominates. 

The remarks we made at the end of Section 13.21 on the 
comparison of the results are also valid for the time delay in 
HS 2209+1914. However, the presence of a possible second so- 
lution in the data and the question of whether the error bars 
should reflect this presence complicate the picture. Again, ba- 
sically diff'erent methods converge towards the same value of 
AtgA ~ 20 days, so a final result of Atg^ = 20 ± 5 is a rea- 
sonable and intermediate solution between possibly too opti- 
mistic and too pessimistic error bars, while not completely ex- 
cluding higher values since these uncertainties are Icr error bars. 
However, given the lack of obvious fast intrinsic quasar variabil- 
ity features in these light curves, and the possible presence of 
microlensing, this time-delay measurement should be handled 
with caution. 



Method 



AtsA in days 



Dispersion 20.99 : 

Regression 20.08 : 

Spline 19.77 : 

NMF 19.28 : 



10.09 (9.74, 2.60) 
13.78(11.50,7.58) 
6.03(5.71, 1.95) 
1.48 



Table 5: Summary of time delays for HS 2209+1914 obtained 
by the four different methods on the complete light curve with 
data from Euler, Mercator, Maidanak, and HCT. The indicated 
error is the total error Between brackets we mention the random 
and the systematic error for the first three methods as explained 
in the text. 



5. Conclusions 

In this paper we present seven-season-long light curves of two 
doubly lensed quasars based on COSMOGRAIL observations 
combining data from Mercator, Maidanak, and HCT telescopes 
for SDSS J1206+4332, as well as with data from Euler for 
HS 2209+1914. All data were reduced in a homogeneous way 
described in [Tewes et al. ( 2012b) (submitted). We applied four 
time-delay estimation methods, based on different principles, on 
these light curves. 

For SDSS J1206+4332 simultaneous deconvolution reveals 
the three galaxies that possibly contribute to the lensing potential 
of the system. We consider the results of the spline fit method as 
our final best time delay of 111 days and pin down the error 
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Subset 






A/g^ in days 
(single run) 


A?g/i± err 
(Monte Carlo approach) 


Mercator 






23.1 


23.20 ± 2.96 


Maidanak 






21.9 


30.48 ±9.32 


HCT 






31.8 


33.32 ± 4.04 


Mercator + HCT 






25.9 


26.50 ±1.95 


Euler + Mercator + 


HCT 




22.7 


26.18 ±1.86 


Euler + Mercator + 


Maidanak + 


HCT 


19.6 


19.28 ±1.48 



Table 6: Summary of time delays for HS 2209+1914 obtained by the NMF on five subsets of data and on the complete set. 

bars to 3.5%, thus improving the uncertainty over a previously 
published result (Paraficz et al.|2009)l. 



Hardly anything has been published on the lensed quasar 
HS 2209+1914 since its discovery. Our work confirms the image 
separation, the flux ratio, and a faint ring-like structure as men- 
tioned in [Chantry et al.| ( |20T0) . We are the first to measure a time 
delay of Af^A ~ 20 days, which was only possible by combining 
data from four telescopes into a 8.5-year-long light curve. 

Even if quadruply lensed quasar systems allow a more 



straightforward estimation of the Hubble Constant Hq, Suyu 
( |2012j ) has shown that doubly lensed systems can augment the 
number of useful time delays when high resolution images of 
the source galaxy in extended arcs can break the radial profile 
slope degeneracy. HS 2209+1914 is promising in this respect 
even if the compactness of the system and the short time delay 
make it difficult to reach sufficiently high precision. 
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